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We study the influence of a high density behaviour of the nuclear 
symmetry energy on a kaon condensation in neutron stars. We find that 
the symmetry energy typical for several realistic nuclear potentials, which 
decreases at high densities, inhibits kaon condensation for weaker kaon- 
nucleon couplings at any density. There exists a threshold coupling above 
which the kaon condensate forms at densities exceeding some critical value. 
This is in contrast to the case of rising symmetry energy, as e.g. for rel- 
ativistic mean field models, when the kaon condensate can form for any 
coupling at a sufficiently high density. Properties of the condensate are 
also different in both cases. 

PACS numbers: 21.65.+f, 97.60.Jd 

1. Introduction 

The possibility that a charged kaon condensate is present in the ground 
state of dense baryon matter has been suggested by Kaplan and Nelson 
[1]. The presence of a kaon condensate would strongly affect astrophysically 
important properties of dense matter in neutron stars [2]. For example, 
the proton abundance of neutron star matter could increase, exceeding the 
direct URCA threshold. This would strongly accelerate cooling of neutron 
stars. Also, the formation of metastable neutron stars could be allowed, as 
the equation of state of hot matter with trapped neutrinos would stiffen, 
supporting larger maximum mass than the cold one [2] . Metastable neutron 
stars with masses exceeding the maximum mass for the cold equation of 
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state would collapse to black holes after the neutron star matter becomes 
transparent to neutrinos. 

The formation of a kaon condensate is expected because of the presence 
of strongly attractive interactions between kaons and nucleons, mainly from 
the so-called sigma term, T> KN . These interactions lower the kaon effec- 
tive mass at higher densities leading to the condensation of kaons above 
a critical density, nk aon , which is estimated in Ref.[l] to be ~ 3no, where 
no = 0.16/m -3 is the nuclear saturation density. 

To study the possibility of a kaon condensation in neutron stars one 
should account for the /3-equilibrium of the neutron star matter which re- 
quires that kaon and electron chemical potentials are equal, fi K - = fi e , 
where the electron chemical potential is given by a difference of proton and 
neutron chemical potentials, /j, e = fip — ^n. This formula shows that also 
nucleon-nucleon interactions, which determine both [ip and /in, are of cru- 
cial importance for the existence of a kaon condensation in neutron stars. 
Thorsson, Prakash and Lattimer [3] have studied the role of nuclear interac- 
tions in a kaon-condensed neutron star matter using simple parametrizations 
of nuclear forces. However, parametrizations used in Ref. [3] are rather 
similar as far as the corresponding nuclear symmetry energy, E sym (nB), 
is concerned. For all models in Ref. [3], E sym {nB) monotonically increases 
with increasing baryon density np. This behaviour at high densities differs 
considerably from that corresponding to several realistic phenomenological 
interactions, such as e.g. \JV\4 + TNI [4], or AVIA + UVII [4], for which 
the symmetry energy actually saturates and then decreases at high densi- 
ties. In our paper we study consequences of this kind of behaviour of the 
nuclear symmetry energy for the formation of a kaon condensate in neutron 
stars. 

The electron chemical potential, fi e , and the proton abundance in the 
neutron star matter are very sensitive to the nuclear symmetry energy at 
high densities [5]. Thus, also the critical density for kaon condensation, 
n kaoni depends on E sym (nB)- The latter quantity is subject to large uncer- 
tainties at higher densities [5]. As indicated above, different models which fit 
the saturation properties of nuclear matter give rather incompatible predic- 
tions of the symmetry energy at high densities. We find that this uncertainty 
affects strongly the critical density of the kaon condensation. In particular, 
for the symmetry energy which decreases at high densities, as e.g. for the 
UV14 + TNI interactions, the kaon condensation is inhibited, at least for 
weaker kaon-nucleon couplings. This is in contrast to the case of monoton- 
ically increasing symmetry energy when the kaon condensate can form for 
any value of the kaon-nucleon coupling at a sufficiently high density. 

In the next section we describe the model of the kaon condensate which 
includes realistic nucleon-nucleon interactions. In Sect. 3 we discuss the high 
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density behaviour of the nuclear symmetry energy. Main results concerning 
the dependence of the condensate properties on the symmetry energy are 
presented in Sect. 4. 

2. The kaon condensate in neutron stars 

The Kaplan and Nelson model [1] of the kaon condensate employs an 
SU (3) x SU (3) chiral lagrangian, obtained in the chiral perturbation the- 
ory, which involves octets of pseudoscalar mesons and baryons. Brown, 
Kubodera and Rho [6] have shown that interactions leading to the kaon 
condensation in the Kaplan and Nelson model are dominated by the s-wave 
kaon-nucleon coupling. Contributions due to remaining interactions are less 
important and we neglect them here. 

When only the s-wave kaon-nucleon interactions are relevant, the kaon 
condensate is spatially uniform and its time dependence is < K~ >= 
vk exp(— ifij^-t) [7], where vk is the amplitude of the mean K~-Reld. For 
such a condensate the effective Kaplan-Nelson lagrangian in the sector in- 
cluding only nucleons and kaons, that is relevant to the neutron star matter, 
reads [3]: 

Lrn = ^u 2 sin 2 9 — 2m 2 K f 2 sin 2 - + n^n(fj,sin 2 - — (2a2 + Aa^)m s sin 2 ^) 

9 9 

+p^p(2fisin 2 - - (2ai + 2a 2 + 4a 3 )m s sm 2 -) + L N . (1) 

Here n = u K - is the kaon chemical potential, 9 = V2vk/ f is the "chiral 
angle", where / = 93MeV is the pion decay constant, and Ljy is the free 
nucleon lagrangian. The parameters, a±, 02, and 03 are coefficients of the 
interaction terms in the original Kaplan-Nelson lagrangian which provide 
splitting of the masses in the baryon octet [1], and m s is the strange quark 
mass. 

The effective lagrangian (1) contains three effective kaon-nucleon cou- 
pling parameters, a\m s , a2m s , and a^,m s (our notation follows that of 
Ref.[3]). The first two are determined by fitting the strange baryon masses 
[3]. In the following we adopt the values aim s = —67MeV and a2m s = 
134MeV from Ref.[3]. The third parameter, a^m s , is related to the kaon- 
nucleon sigma term, 

S XAr = --(ai+2a 2 +4o 3 )m s , (2) 

which is poorly known and thus the value of a^m s is subject to a consid- 
erable uncertainty. We use here values in the range —134MeV > a^m s > 
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—3WMeV which correspond to the sigma term in the range 170MeV < 
Y, KN < 520MeV [3]. 

The lagrangian (1) describes only the kaon-nucleon interactions in the 
condensate. To account for nucleon-nucleon interactions, which are crucial 
in a dense nucleon matter in neutron stars, we use realistic models of the 
nucleon matter [4] . The energy density of the neutron star matter with the 
kaon condensate is composed of three contributions, 

e ns = £kn + CN + Qep, (3) 

where 6kn is the energy of the kaon condensate described by the lagrangian 
(1), e^v includes both the nucleon-nucleon interaction contribution and the 
nucleon Fermi energy, and e; e „ is the electron and muon contribution. 

The energy per particle, ew/ng, obtained in variational many-body 
calculations with realistic nucleon-nucleon interactions, can be written as 
a function of the baryon number density, ng, and the proton fraction, 
x = np/riB, in the form [8] 

E(n B ,x) = T F {n B ,x) + V {n B ) + (1 - 2x) 2 V 2 (n B ), (4) 

where Tp(n B ,x) is the Fermi-gas energy, and Vo(n B ) and ^("b) are the 
interaction energy contributions. 

The energy density of the neutron star matter with the kaon condensate 
including nucleon-nucleon interactions reads 

e ns (n B ,x, fi, 9) = T F (n B ,x)n B + n B V (n B ) + n B (l - 2x) 2 V 2 (n B ) 

9 fi 2 9 
+2m 2 K f 2 sin 2 - + /— sin 2 9 + (2aim s x + 2a 2 m s + 4a^m s )n B sin 2 - + e e + e M , 

(5) 

where the electron Fermi sea contribution is 

£e = h- (6) 

Muons are present only when the electron chemical potential, /x, exceeds 
the muon rest mass, fi > m^. The energy density of the muon Fermi sea is 

= m^fiPFjm^), (7) 

with 



The above formulae allow us to determine the critical density of kaon 
condensation and to study properties of the neutron star matter with a 
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developed kaon condensate. We obtain the ground state parameters, the 
proton fraction, x, the kaon chemical potential, /j,, the kaon condensate 
amplitude, 9, and the energy density, e ns , as functions of the baryon number 
density, ng, by optimization of the thermodynamical potential e = e ns — 
l J '(. n K + n e + n M ~~ n p) which reads 

e(n B ,x, /i, 9) = T F (n B ,x)n B + n B Vo(n B ) + n B (l - 2x) 2 V2(n B ) 

i-2^ 2 ■ 2 n , o 2 f 2 ■ 2 ^ , /i , \ • 2 ^ 

— / — sm 9 + 2m K j sin - + fj,n B x — /me(i + xj sin - 

+ (2aim s x + 2a 2 m s + 4a 3 m s )n B sin 2 - + e e + e M . (9) 

This thermodynamical potential is related to the (Landau and Lifshitz) 
potential, 17/ V = e — S/ijnj = e — fj>Nn B , where /Uat is the neutron chemical 
potential. 

The critical density, nfc aora , is defined as a density at which the con- 
densate amplitude starts to deviate from zero. Our procedure provides the 
chiral angle of the condensate as a function of density, 9(n B ). The value of 
rikaon is thus found from the condition 9(rikaon) = 0. 



3. The nuclear symmetry energy at high densities 

The nucleon-nucleon interaction energy is parametrized in Eq.(4) in 
terms of the isoscalar and isovector contributions, Vo(n B ) and V^n^). As 
one can notice, at a given baryon density, n B , the interaction energy den- 
sity corresponding to the isoscalar part, n B Vo(n B ) in Eq.(9), is constant and 
thus it does not play any role in the optimization of the potential e which 
determines the kaon condensate parameters. It is the isovector contribution, 
V2(n B ), which is crucial for the onset of the condensation. This component 
is directly related to the nuclear symmetry energy which expressed in terms 
of V2(n B ) reads 

E sym (n B ) = —Tf{ti b , -) + V 2 (n B ). (10) 

As we already mentioned, the high density behaviour of E sym {n B ) is not 
well known at present. Different model calculations give incompatible ex- 
trapolations away from the empirically determined value at no, E sym (no) = 
31 ± AMeV [9]. In some approaches [5], E sym (n B ) is a monotonically in- 
creasing function of the baryon number density, whereas other models [5] 
predict E sym (n B ) to saturate and then decrease at high densities. Physi- 
cally, in the former case the energy of pure neutron matter is always higher 
than the energy of symmetric nuclear matter, while in the latter one pure 
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neutron matter becomes eventually the ground state of dense baryon mat- 
ter. This difference has profound consequences for neutron star matter. 
Here we study how it influences the onset of the kaon condensation. 




2 4 6 8 10 



n/n 

Fig. 1. The nuclear symmetry energy as a function of density for indicated models 
of nuclear interactions. The curve E\ corresponds to a linear parametrization of 
Ref.[3]. 



In Fig.l we show the symmetry energy corresponding to three different 
models of nucleon matter which we use in our calculations. The curve 
labelled E\ corresponds to a linear parametrization used in Ref. [3] . It fits the 
empirical value, E sym (no), and then extrapolates linearly to higher densities. 
Such a linear dependence of the symmetry energy on density is typical 
for relativistic mean field models of nucleon matter [5]. Other curves in 
Fig. I show the symmetry energy from Ref. [4] obtained in variational many- 
body calculations with realistic two-body nucleon-nucleon potentials and 
three-body interactions. The curves labelled UV1A + TNI and AV1A + 
UVII correspond to the Urbana v 14 two-body potential with the TNI three- 
body term [4] and the Argonne i>i4 potential with the UVII three-body 
interaction [4], respectively. As one can notice, realistic models predict that 
the symmetry energy saturates and then decreases with increasing baryon 
density. This behaviour is markedly different from that of the curve E\. 
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4. The critical density and properties of the kaon-condensed 



In this section we present results for three models of the nuclear sym- 
metry energy shown in Fig.l. Particularly interesting is the corresponding 
critical density, rikaon- To best illustrate the influence of the symmetry en- 
ergy on rikaon we use a different definition of the critical density which states 
that rikaon is the lowest density for which the energy of the lowest state of 
K~ in the neutron star matter, w_, becomes equal to the electron chemical 
potential, 



At higher densities, n > rikaon, kaons form a Bose-Einstein condensate of 
finite amplitude. This definition of the critical density is equivalent to that 
given at the end of Sect. 2. 

The lowest kaon energy in the nucleon medium reads: 



The energy to- decreases with baryon number density, Figs. 2 and 3. The 
condition (11) can be satisfied only at a sufficiently high density provided the 
electron chemical potential does not decrease too fast. In Fig.2 we show the 
lowest kaon energy, w_, for three values of the coupling parameter a^m s , 
and for the electron chemical potential, [/,, corresponding to the UV1A + 
TNI interactions. The electron chemical potential is rather low, it has a 
maximum at n = 3no and then decreases to zero at a density of about 6wo, 
where electron and proton densities vanish. The critical condition (11) can 
be satisfied only for the coupling values at least as strong as in case (b), 
a-zrn s < — 230MeV. For weaker coupling the u>- curve does not cross the 
chemical potential and the critical condition (11) can never be satisfied. 
Hence, in the case of the symmetry energy of the UV1A + TNI potential 
there exists a threshold value of the kaon-nucleon coupling parameter, a^m s , 
below which there exists no kaon condensation in the neutron star matter. 

Such a threshold does not exist in Fig. 3 where we show the kaon energy 
cj_ and the electron chemical potential for the symmetry energy E\. Since 
the electron chemical potential monotonically increases with density, the 
critical condition (11) can be satisfied for any value of the coupling param- 
eter a^rris at a sufficiently high density. Critical densities corresponding to 
indicated values of a^m s in Fig.2 differ by a factor of about two. 

Values of the critical density can also be read off from Fig.4 and Fig.5, 
where we show the amplitude of the condensate, 0, as a function of density, 
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Fig. 2. The lowest kaon energy in the nucleon matter, (solid curves) and 
the electron chemical potential, /x, (dashed curve) for the UV14+TNI interac- 
tions. Curves labelled a,b, and c correspond to the coupling parameter a 3 m s of 
-ZIOMeV, -230MeV, and -134MeV, respectively. 
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Fig. 3. The same as in Fig. 2, for the £1 symmetry energy 



for the UV14: + TiVT and AV14 + UVII interactions, respectively. In both 
figures, results corresponding to the symmetry energy E\ are also shown. 
The values of rikaon found from Fig.2 and Fig.3 coincide with those from 
Fig.4. It is also interesting to compare how the angle 9 varies with density. 
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In Fig.4, for the UV14 + TNI model the amplitude grows very fast with 
density for densities above the critical one. It reaches a maximum when 
the proton fraction becomes x = 1, and decreases at higher densities. A 
similar behaviour is displayed in Fig. 5 for the AV14 + UVII interactions 
for the strongest coupling (solid curve a) . This is in contrast to the E\ case, 
when the condensate amplitude monotonically grows with density for all 
indicated values of the coupling parameters. 




2 4 6 8 10 

n/n 

Fig. 4. The condensate amplitude as a function of density for the UV1A + TNI 
interactions, solid curves, and for the symmetry energy Ei, dashed curves. The 
coupling parameters corresponding to labels a and b are the same as in Fig. 2 

In Fig. 6 the proton fraction of the neutron star matter with the kaon con- 
densate is shown. One can notice that for the UV14 + TNI interactions the 
kaon-condensed phase contains only protons. Neutrons are fully converted 
into protons at densities only slightly exceeding the critical value. This 
behaviour is similar for all coupling parameters stronger than the thresh- 
old value 03777.3 = —230MeV . In Fig. 7, where the same is shown for the 
AV14 + UVII interaction, one can notice that for the strongest coupling 
the proton fraction also reaches unity at high densities. For interactions 
with the symmetry energy E\ the proton fraction at high densities tends to 
an asymptotic value x = 0.6, for all kaon-nucleon couplings. 

5. Discussion 

The above results show that the critical density of kaon condensation in 
neutron stars is rather sensitive to the high density behaviour of the nuclear 
symmetry energy. For some realistic equations of state, such as the UV14 + 
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Fig. 5. The same as in Fig. 4 for the AV14 + UVII interactions. 




n/n 

Fig. 6. The proton fraction of the neutron star matter with a kaon condensate for 
the UV\A + TNI interactions and for the symmetry energy E±. The curves a and 
b correspond to values of the coupling parameter given in Fig. 2 



TNI one, the existence of the condensate is even forbidden if the kaon- 
nuclcon interaction is not strong enough. Also, properties of the neutron 
star matter with developed condensate are sensitive to the symmetry energy. 
One can notice in Fig.4 and Fig.5 that the condensate amplitude, measured 
by the angle 6, behaves in a quite different way for considered models of 
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nucleon-nucleon interactions. In particular, the kaon condensate weakens 
at high densities for the UVIA + TNI interactions, Fig.4. 

For astrophysical applications, the proton fraction of the neutron star 
matter plays an important role. The influence of the symmetry energy 
on the proton abundance in the kaon-condensed phase is rather strong, as 
shown in Fig.6 and Fig.7. In fact, the proton abundance at higher densities 
is determined entirely by the nuclear symmetry energy. The kaon-nucleon 
coupling parameter, a^m s , affects the proton fraction only at densities close 
to the critical value. The role of the symmetry energy is most spectacular 
for the UV14 + TNI interactions for which neutrons are fully converted into 
protons, Fig.6. 
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Fig. 7. The same as in Fig.4 for the AVU + UVII interactions. The curves a and 
c correspond to values of the coupling parameter given in Fig. 2 

Let us remark finally that recent claims [10] that all modern phase equiv- 
alent potentials which fit accurately the n — n and n—p scattering data yield 
a symmetry energy which increases with density, are unjustified. For the 
Argonne potential AF18 authors of Ref.[10] find increasing E sym (nB)- Cal- 
culations reported in Ref.[ll] show that the symmetry energy corresponding 
to this potential saturates and then decreases at high densities in a similar 
way as found by Wiringa, Fiks and Fabrocini [4] for Urbana UV14 and 
Argonne AV14 potentials. Premature conclusions of Ref.[10] stem from the 
use of the lowest order Brueckner approach which is inadequate at high 
densities. Hence the uncertainty as to the high density behaviour of the 
nuclear symmetry energy is still present. For astrophysical applications, 
both decreasing and increasing symmetry energy should be used in order to 
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assess the role of this uncertainty. 
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